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ABSTRACT 

The diffraction by a material discontinuity in a thick dielectric/ferrite layer 
is considered by modelling the layer as a distributed current sheet obey- 
ing generalized sheet transition conditions (GSTCs). The sheet currents 
are then formulated and solved via the standard dual integral equation ap- 
proach. This yields the diffracted field in terms of unknown constants which 
underscore the non-uniqueness of the GSTC current sheet representation. 
The constants are dependent on the geometry and properties of the discon- 
tinuity and are determined by enforcing field continuity across the material 
junction. This recpiires the field internal to the slab which are determined 
from the external ones via analytic continuity. Results are given which 
validate the solution and demonstrate the importance of the constants. 


OBJECTIVE 


This task involves the use of higher order boundary conditions to generate 
new solutions in diffraction theory. In particular, diffraction coefficients will 
be developed for dielectric/magnetic layers and metal-dielectric junctions 
which are often encountered on airborne vehicles as terminations of coatings 
and conformal antennas. Solutions for both polarizations will be developed 
for fairly thick junctions and versatile computer codes will be written and 
tested. Creeping wave diffraction coefficients will be also developed for 
multilayered coated cylinders. 


PROGRESS 
1 Introduction 

In scattering, layered materials are often modeled by equivalent sheets 
satisfying simple boundary/transition conditions. In particular, impene- 
trable layers are typically replaced by opaque sheets satisfying standard 
impedance boundary conditions (SIBCs) [1], whereas penetrable layers are 


13 


PRECEDING PAGE BLANK NGT FILMED 



represented by transparent sheets obeying resistive or conductive type tran- 

at" rno' 0n Vfi U / P i TheSe SimP ' e b “ a <>ary/transition conditions 
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solution constants [12] and this simply points to the fact that additional 

conditions are required for their specification. 

In this paper we demonstrate that the GIBC/GSTC sheet characteriza- 
tion can yield a complete solution when supplemented with certain condi- 
tions at the sheet discontinuity which do not require apriori knowledge of 
the edge fields. As a vehicle in presenting this solution procedure we employ 
the dual integral equation method to consider the plane wave diffraction 
by a discontinuous distributed sheet (see Figure 1(b)). This very general 
model is capable of representing material half-planes, material junctions, 
and material discontinuities on grounded structures, such as those shown 
in Figure 2. In addition, a distributed sheet model typically renders the 
same degree of accuracy as the usual infinitely-thin sheet, but with a lower 
order condition. It is, therefore, of much practical interest. 

In the first part of the paper, the GSTC representation of the distributed 
sheet discontinuity is used to develop dual integral equations in terms of 
the unknown spectral functions proportional to the sheet currents. These 
equations are then solved in the standard manner to yield expressions for 
the spectral functions in terms of unknown constants, and examples are 
presented where a proper choice for the constants demonstrates that they 
recover known solutions. This demonstrates the validity of the presented 
solution, but in general, the determination of the constants requires the 
enforcement of additional constraints demanding field continuity across the 
junction. The development of these conditions and their use in solving for 
the constants is also presented. 


2 Dual Integral Equation Formulation 

Consider a distributed sheet of thickness t illuminated by the plane wave 


F — c 

1 inc — ° 


jk(x cos <t>o+y sin <j> 0 ) 


E Z ' inc , E z polarization, 
ZoHz'Xnc H* polarization, 


( 1 ) 


as shown in Figure 1(a). The excitation (1) induces reflected and trans- 
mitted fields which are explicitly given by the properties of the distributed 
sheet. If this sheet models a symmetric slab, then an appropriate GSTC 


15 



representation is formally given by [10] 

u '" (^) _ • F "} + k U “ (~fr) [*" " + *“]} = 0, 

— oo < X < OO 


«« I - 


&){** + *-} + {dy [r - f -] } = o, 


-oo < x < oo. 


( 2 ) 

in which F is the total field, F* = F(x,y = ±r/2). = Ah, „ 

_ ±T 2) ’ a u d ^ F± = Also, ’4 (_£) are differential 

finite ooT r rat l°" the SeW quantily ” lhe brackets, and are 

te polynomials in whose coefficients depend on the slab modeled 

by the distributed sheet. To maintain the generality of the solution, the lit. 

p xa ors are left in symbolic form and the reader is referred to [10] for their 

_ p cit representation in terms of the material constants and thickness of 

the ZT ™ m?r T S the m ° dded Slab ‘ In gCneraI ’ the order of U\ x (i e 
the highest derivative present) is usually the same or one more than that 

off U and similarly the order of U\ x is the same or one more than the order 

oi U 22 . Thus, we may define the orders of the GSTCs in (2) to be 

N ° dd = maximum { order of U\ x (A 2 ) , 1 + order of U\ 2 (A 2 )} 


jyeven 


f (3) 

- maximum { order of U\ x (A 2 ) , 1 + order of U\ 2 (A 2 ) } 
emp^" fi“d ‘ IanSmitted fieWS ^ n ° W be eaS “ y ^ 


Frefl = R x e 3k(xcos<t> °-y sm <t>o) 
Ftran = T x «**<>+» sin *„) 


(4) 

(5) 


r d fZ the reflection “ d “ sio " 


gjArr sin <f> Q 

R x = + R odd] 


T\ = 


2 

gjArr sin <j > 0 
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[R e r n - R° x dd \ . 


( 6 ) 

(7) 
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with 


sin <j>M\ 2 ( CQs2 $ 0 ) - U\\ ( cos2 4*°) 
sin 4>M22 (cos 2 4 >o) + ( cos2 0 °) 

sin ^Mi 2 (cos 2 <fto) ~ ^11 (cos 2 <ft Q ) 
sin <j>JU\2 (cos 2 <f> 0 ) + U\ x (cos 2 <f>„) 


(S) 

( 9 ) 


We remark that in (8) and (9), U} } (cos 2 <p 0 ) now represent simple polynomial 
functions in cos 2 <j> 0 , since -dx 2 /k 2 = cos 2 <f) 0 in view of the field expressions 
(4) and (5). 

Consider now the case where the right half of the distributed sheet 
in Figure 1(a) is replaced by another sheet of the same thickness, but of 
different properties, as illustrated in Figure 1(b). The GSTC representation 
of this modified sheet is 


( 10 ) 

for —00 < x < 0 and 

(11) 

for 0 < x < 00 , where the superscripts 1 and 2 distinguish the left- and 
right-hand sheets, respectively. Referring to our previous discussion, the 
orders of the right hand side GSTCs are given as 

J\rodd — max j order of U xl in A, 1 + order of U x2 (^ 2 ) m 

(12) 

N' ven = max {order of U\ x (A 2 ) in A, 1 + order of U\ 2 (A 2 ) in a} . 
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The modified right hand side sheet induces a scattered field F s in the 
piesence of the excitation (1), and the total field can be represented as 

p _ / Fine + F re fi + F 3 y > t/ 2 

\ Ftran + F, y < r/2 ' ' 

where F s is the unknown scattered field in the region |y| > t j 2 and can be 
expressed as [13] [14] 


F'(x,y) = f 

Jc 


\y\ 


Podd (cos a) + P even (cos a) 


e -jksina(\y\-T/2) e -jkxcoso, ( f a ^ -j^ 

where C is a contour in the complex a plane, such that A = cos a runs from 
— oo to oo as shown in Figure 3. In this, the spectral functions P 0 dd (cost*) 
and P even (cos o) are directly related to the Fourier transforms of the un- 
known equivalent currents 


Jodd = F s + - F s (15) 

Jeven = F+ + F g , ( 16 ) 

via the relations 

/ °° d \ 

Podd (A) e~ jkxX —r=== (17) 

-CO y/l — A 2 

/ oo JX 

P even (\)e-^ X ~ A (18) 

-oo v 1 ~ A 2 

Substituting (1), (4), (5), (13) and (14) into the transition conditions (11) 
and (12), and introducing the transformation A = cost* (see Figure 3) yields 

/“er(V)p^( A)e-^-^= = 0 (19) 

|“ ar ,( A 2)p t<jm(A)e - illJ _£_ = 0i (20) 
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for x < 0 and 


2sm<f> 0 e lkxX °e jkT/2s,n *°Z odd (\l) 


Qr( A») 

_ 2 sin sin 6 °Z ei;e7l ( ) 

^eve„( A 2) 

for x > 0, where A 0 = cos<p 0 and 

Q\ dd { A 2 ) = u l n (a 2 ) + v/T - A m\ 2 (a 2 ) 

geven ^ = ^ ^ + ^ _ A 2^ ^2^ 

SI" (a 2 ) = w?, (a 2 ) + >/i - x 2 u 2 X2 (a 2 ) 

$r n (a 2 ) = W|, (a 2 ) + n/ 1 - A 2 z4 (a 2 ) 


( 21 ) 


221 


(23) 

(24) 

(25) 

(26) 


Zm(K) = 

U'n (Aj)«?a (a*) -UU ($«?■ (Aj)J 

(27) 

Zeven (A^) = 

K(Aj)«l,(Aj)-«i(A;)M|,(^)]- 

(28) 


Equations (19) with (21) and (20) with (22) form two uncoupled sets 
of integral equations, sufficient to yield a solution for the unknown spectra 
Podd(^) an d P even (X). Clearly, because of the similarity between the two 
sets of equations, once a solution for P od d(X) is found, the corresponding 
one for P even (A) follows by inspection. 
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3 Solution of the Diffracted Field 


Upon a solution of the dual integral equations (19) and (21) we obtain 

_ \2 e jfcr/2sin</) 0 


j sin <j>, 

PmW ~ 


A n 


godj (\)godd {K) godd (A) godd (Ao) 



Eodd (A) 

E 'odd ( A n ) 


N 0 dd. — 1 Nodd~ 1— TO 

+ XI XI °mn (A + A 0 ) m (AA 0 ) n 

m— 1 n=0 


(29) 


where_we have assumed that J 0 dd(x) ~ x s ° dd as x — ► 0 with 0 < S od d < 1- In 
this, N odd = int {1/2 {Nl dd + Nl ld + 1)}, and a mn are arbitrary constants 
as yet undetermined, and correspond to the coefficients of the polynomial 
resulting from the application of Liouville’s theorem. The chosen symmetric 
form of this polynomial is not unique but will be found most useful later in 
constructing a reciprocal form for P od d{ A). Also, G 1+ (A) are Wiener-Hopf 
split function regular in the upper ( + ) or lower (-) half of the A-plane and 
satisfy the relation 


6h(A 2 ) = G 1+ (A)Gh_(A) 


(30) 


(see Appendix). Similarly, C? 2 ±(A) are the corresponding split functions 
associated with (^(A 2 ). Finally, E odd ( A) is some entire function behaving 
no worse than |A|(^ dd+jV °^) /2-Sodd and can take any of the forms 


Eodd (A) = 


" Z Q dd (— AA 0 ) 

ztddW 


or 

or 


(31) 


where Z ± ( A) are again upper and lower functions satisfying the relation 


Z(A 2 ) = Z + (X)Z~(\) (32) 

Following a similar procedure we obtain P even (A) as 
P /> ) _ j sin (f> 0 y/l - A 2 e jkT/2sm<t>o 

even - 2?r a + a 0 (A (a 0 ) (A) g fip ( a 0 ) 


( \ 2^ ^ even 

*°) F 


/ ^ A ^et«n 1 Neven — 1 — TO 

JZT-+ E E ^(A + AJ" , (AA„)“ 

\ A o) m = i 


71 — 0 


(33) 
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with E even (A), N even and b mn being the counterparts of £w^(A), N 0 dd and 
a mn , respectively. Taking into account the choices (31), we may substitute 
(29) and (33) into (14) and subsequently perform a steepest descent path 
evaluation to obtain for p — » oo (all surface wave contributions are neglected 
in this evaluation) 


e ~3 k P 

F{p,<i>) ~ [D 0 dd{<i>,<t>o) + D even {4>,<f> 0 )\ . (34) 

where (p, 0) are the usual cylindrical coordinates and D a dd (0, (j> 0 )+D e V€n (0, <f> 0 ) 
is the far zone diffraction coefficient symmetric with respect to 0 and 0 O . 
We have 


D Q dd (0? 00.) — 




sin 0 o sin 0 


27T cos 0 4* cos 0 O 

gjfcr /2(sin<£ 0 4*| sin 0|) 


Deven (0* 0o) — 


(cos 0) fffr (cos 0 O ) ( COS 0) ^ 2 + ( COS 0o) 

^ odd 1 ^ 1 ^ 

(cos^,cos<i» 0 ) 4- E E ^mn 

m = l n=0 


(COS 0 +COS 0 o ) m (COS 0 COS 0 o ) n 

e” J7r / 4 sin 0 O | sin 0| 

2tt cos 0 + cos 0 O 

gjkr/ 2(sin^ 0 +| sin<^|) 


(35) 


^i e H en (cos 0) £/f!. en (cos <^ 0 ) Q^ n (cos <f>) Q e ^ n (cos <f> 0 ) 

^ ^even 1 Neven 1 171 

^ even (COS 0, COS 0o) ^ ^ ^ ^ ^mn 

m=l n=0 


• (COS 0 +COS 0 o ) m (COS 0 COS 0 o ) n 


(36) 
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ill which the functions Z 0 dd,even (cos <f>, cos <f> 0 ) are given by (see (31)) 

{ z odd ( - COS <f> cos <t>o ) or 

Z odd( cos< t>) Z odd( cos< f>o) or (37) 

Z odd (cos 4>) Z+u (cos <f>o) 
f Z even (— cos <j> cos <j>o) or 

Z even (COS (j) , COS 4>o) = { Z tx ,en ( COS &) Z even ( COS ^ o) OT (38) 

. Z even ( COS <f>) Z + en (cOS (j> 0 ) 

Because the above three choices for Z 0 dd and Z even differ only by terms of 
the form (cos (f> + cos <j) 0 ) m (cos <f> cos 4 > 0 ) n , it is immaterial which of them we 
choose, although one of the choices may likely lead to a more compact rep- 
resentation. Nevertheless, regardless of the choice of Z a dd and Z even , one is 
still faced with the determination of the unknown constants a mn and b mn in 
(35) and (36), repectively. These are a manifestation of the non- uniqueness 
of the finite-order GSTC sheet model employed herein, and their explicit 
determination requires the introduction of additional constraints pertain- 
ing to the physics of the problem. Before we consider their determination 
for the general case, we first look at a specific example, that of diffraction 
by a thin single layer junction. 

4 Diffraction by Thin Single Layer Discon- 
tinuous Slabs 

The diffraction coefficient given by (35) and (36) is very general and can 
model a wide variety of geometries. To check its validity, display its ver- 
satility, and assess the relative importance of the unknown constants, we 
consider the thin material-to-material junction of thickness 2 w as shown 
in Figure 4. The slab will be modelled by a sheet of thickness 2 (w — w 3 ) 
and with a proper choice of the material parameters this geometry can re- 
duce to junctions whose diffracted field is available, thus, permitting some 
validation of our solution. 

If the left hand side of the slab, in addition to being thin, is also asso- 
ciated with low index of refraction, it may be modeled by a low contrast 
GSTC sheet. Thus, an 0(ud,iu]) approximation with terms of 0(w 3 w ) ne- 
glected is sufficient for the representation of the operators or polynomials 


22 



U\ . In particular, we have 


U\ x (- dx 2 /k 2 ) = 1 

u\ 2 (— (9i 2 /t 2 ) = j k. (u\W — w s ) 

u\\ {-dx 2 jk 2 ^ = jk ( ** - u; 5 ) + ^ w^jdx 7 

U\ 2 (~dx 2 /k 2 ) = 1 


( 39 ) 


where ej and f(\ are the relative permittivity and permeability of the left 
hand slab, respectively, and 


ill 


fi i, E z polarization 
€i H z polarization 


(40) 


Also, when w s = u% these are simply the transition conditions derived first 
by Weinstein [5] and later by Senior and Volakis [6]. The corresponding 
polynomials to be employed in (23) - (28) are given by 

Uh ( — cos <f) cos 4> 0 ) — 1 

U\ 2 ( — cos <t> cos (p 0 ) =z jk ( lliW — w s ) 

U\ x (- cos cos <f> 0 ) — jk ( J) 1 - u> 3 ) + jk u> 3 ) 

lll 2 ( — COS (j) COS <j>o) = 1 


COS 4> COS <po 


Incorporating these into (35) and (36) and setting 

Z 0< id (COS (p, COS <f> 0 ) = Z a dd ( — cos <p COS (f> 0 ) 


Zeven (COS <^, COS <f> 0 ) = Z even ( - COS <f) COS <f> Q ) 


(41) 


(42) 

(43) 


yields 


= - e ~ s i n sin ^ 

27 r cos cp + cos (p 0 

1^12 ( — COS <f) COS 4> 0 ) — jk(u I W — IV s ) U xl ( — cos <f> cos <f> 0 ) 
A/_ (cos <j>; Y dd,x ) M- (cos <p 0 ; 7 odd>1 ) §%+ ( cos 4>) £%+* ( cos 4>o) 


(44) 
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D even (<^ 5 & 0 ) 


e~ J77 / 4 sin<ft 0 | sin<ft| jfer/2(sin 4 > 0 -H sin<£|) 
27T COS 0-f COS <i>o 


a 3 


[o 1 COs 0 cos M 2 I — COS 4> COS <fto) — ^21 ( — COS <f) COS 4>o ) 

[rim-1 M-( y COS( t > i'Yrn )M- (cos i^o^m )]^2+ Cn ( COS ^)^2+ Cn ( COS ^°) 


610 (cos (ft+COS (ftp) 


w 


+ ajnl =1 M-(costa;""'' 1 }W-(c«?°;im' n '' )]C|;-(cos0)S? 2 --(cos0 o ) 
/here the split fu "'* 1 'on M_ (cos <f>; 7 ) is given in the Appendix, 


-7 (45) 


and 


odd, 1 


et'en,l 

7l,2 


a x 

(we ifii 

- jk I w. 

\ U\ 

<*2 

= ik (vr w ) 

Oi 3 

= jkW {e m -l) 

U 1 


-3 

k (u 1 

w - w s ) 

U 1 ± 

\]u\ + Ak 2 w(e l f.ii - 


2j A: (u/ — 


(46) 


(47) 


with 7 0<i(i or eue " are associated with possible surface wave poles. To complete 
the definition of ( 44 ) and ( 45 ), the functions associated with the right hand 
side properties of the slab (i.e. those functions with the superscript 2 ) must 
be specified and Tables 1 and 2 provide explicit expressions for the functions 
Ufj (-cos^cos <j> 0 ), Q°£ (cos <f>) Q° 2 d + d (cos 4> 0 ) and Q^ en (cos </>) Q e 2 v + en (cos (j> 0 ) 
terms. By edge condition considerations, all of the constants Qmn and b m n 
have been set to zero except 610 appearing in the definition of Devem which 
is non-zero unless the right hand side slab is a PEC/PMC under an E z j H z 
excitation (see Table 2). 

By invoking image theory the diffraction coefficient for the 
grounded metal- dielectric join, shown in Figure 5 is given by 

E^ rs (cOS <j), COS (j ) 0 ) = 2 D eue n(cOS <f>, COS (j> 0 ) 
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The GSTC or GIBC model for this structure cannot dsicriminate whether 
the stub at the junction is a perfect electric conductor (PEC) or perfect 
magnetic conductor (PMC). This information can only be carried by the 
constant b\o and its determination must somehow involve the properties 
of the junction across its thickness as discussed in a subsequent section. 
However, since the diffraction coefficient for the junction in Figure 5 is 
already available [15], &io can be identified. Upon setting w s = 0, we find 


^pec stub 




r-i)[M-(> 


r - ' ewcn.l \ I # / / * ■ * ■ c v 

/ u7uo 1 )M-\V C i^io 2 



j kw 


0( A '-C 


iiL 
U 

ev«n,l 


)m. 


. even 

/€ 1 V 1 0 2 




(49) 

(50) 

(51) 


This comparison clearly demonstrates the importance of the constant &io 
and by referring to Figure 6 we observe that it plays a major role in the 
computation of the diffracted field. 


5 Modal Decomposition of the Symmetric 
Slab Fields 

A general approach for determining the solution constants is to enforce 
tangential field continuity accross the junction. This, of course, demands 
a knowledge of the fields internal to the discontinuous slab, which are not 
readily available when a GSTC simulation is employed. The Weiner-Hopf 
(or dual integral equations) solution in conjunction with the GSTC pro- 
vides only the fields external to the slab, and this section deals with the 
determination of the internal fields from the external ones. 

A modal representation of the internal field is first proposed comprised 
of discrete and continuous spectral components. This representation is 
compatible with that given by Shevchenko [16] whose eigenfunctions are 
chosen to satisfy field continuity across all layer interfaces including the 
air-dielectric interface. Consequently, the representation is valid inside and 
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outside the dielectric once the coefficients of the modal representation are 
determined. This is accomplished by recasting the Weiner-Hopf or dual 
integral equation solution given earlier (see (13), (14), (29) and (33)) in a 
form compatible with the proposed modal representation, thus permitting 
the identification of the modal or eigenfunction coefficients. These will, of 
course, be in terms of the unknown constants appearing in the Weiner-Hopf 
solution and the enforcement of field continuity accross the junction leads 
to a linear system of equations to be solved for the constants as described 
in the next section. 

For the symmetric slab in Figure 2, the total field may be decomposed 
into its odd and even components. Specifically we write 


F (*,y) 


F Uodd (x,y) + F 1 ' euen (x,y) x<0 


F 2 ' odd (x,y) + F 2 - even (x,y) x>0 

where F odd (x,y) = — F odd (x, — y) and F even (x,y) = F even (x, — y). Follow- 
ing [16], the odd and even fields interior and exterior to the slab may be 
expanded into discrete and continuous eigenmodes as 


,odd 


{x,y) = 


m = 1 7 




+ E B]^ >dd ^]^ dd (y ) e ~ jkxX ™° dd 

m= 1 

roo . . 

+ J C hodd (f3) V 1 ' odd (A 2 , y) e~ ]kxX d0 (53) 


'll, even 




,-jkx \'J° 


90 / 2 \ 

E 4,r n *i'' v “ {(w ) ,») 

m= 1 ' ' 

A; 1 . even 

* V 3U> 

+ E B l m ewn ^ m even (y)e-^ 

m = 1 

+ jT° (/?) ^ lieven (a 2 , y) e~ ikxX dp 


\ 1 .even 


( 54 ) 


F 2,0<w (x, y) 


E-4^^((A^f,y) 

m= 1 V 7 


e 


—jkx X 


2, go 
m 


26 



Y 2, odd 

ly 3W 


| ^ j^2,odd q 2, odd 


(y) e 


—jkx A 


2, odd 
m 


m — t 


roc , x 

+ j C 2 ' 0<w (0) V 2 ' odd (a 2 , y) e~ jkxX d0 


( 55 ) 


^2,ei/en 





Y2,<ven 

+ £ 

771= 1 

/“CO / x 

+ j C— (0) $ 2 ' even (A 2 , y ) e~ ]kxX d0 


B l 


(y)e-> kxX 




(56) 


where /m jA^’ 6 " 6 "} < 0 and A = \/l — /? 2 , with the branch of the square 
root chosen so that Im{\J 1 — /? 2 } < 0. In (53)- (56), <I> odti ' evcn are referred 
to as the cross section functions corresponding to the continuous modal 
fields whereas <$° dd ' even are the corresponding cross section functions for the 
discrete modal fields associated with the surface waves. The cross section 
function associated with the geometrical optics fields is also \^ odd ^ V€n eval- 
uated at A = A^°, where A^°is a parameter to be determined later. As can 
be observed from (53) - (56), the cross section functions specify the field 
behavior in the plane normal to the slab, and hence all information per- 
taining to the fields interior to the slab are embedded into these functions. 
They will be chosen to satisfy the orthogonality relations (where u(y) is 
l-i(y) or e(y) for E z or H z polarization, respectively) 



fr(A 2 ,y)*(A 2 ,y) 

«(y) 


dy = 0 for A ^ A 



$ m (y)*(A 2 ,y) 

«(y) 


dy 


0 


(57) 

(58) 


and thus each discrete eigenmode <3> m (y) e~ jkxX and each continuous eigen- 
mode ( A 2 , y ) e~i kxX must satisfy the wave equation. Additional details 
pertaining to the cross section functions are given in [16]. 
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Exterior Cross Section Functions 


To compute the cross section functions in the exterior slab region |y| > 
r/2, we recall that in accordance with the slab simulation based on the 
generalized sheet transition conditions (GSTCs), the external fields satisfy 
the conditions (10) and (11). Because of the orthogonality relations (57) 
and (58), each of the cross section functions (A 2 , y) and (y) must then 
satisfy their respective odd or even GSTC. In view of this we set 



T7/j2 (a 2 ) cos [k(\y\ - r/2) n/T^A 2 ] j 

= Y {^n (A 2 ) <7 i 2 (1,T M -r/2, A 2 ) 

+U V X2 (A 2 ) <722(1, 1, \y\ - r/2, A 2 )} (59) 



+^22 (^ 2 ) cos [ k (|y| - r/2) Cl - A 2 ] } 

= {^21 (^ 2 ) 912(1, 1, |y| - r/2, A 2 ) 

+U 22 (A 2 ) 922(1, 1, |y| - r/2, A 2 )} (60) 

where yp represents the infinite order form of the q tJ layer operators given 
in [10]. Once each of the modes comprising (53) through (56) is substituted 
into (10) or (11), the differentiation implied by —dx 2 /k 2 reduces to a mul- 
tiplication by A 2 and the above 'I '?’ odd and <&P’ even are then readily shown 
to satisfy the associated GSTC. It can also be shown that these satisfy the 
orthogonality conditions (57) and (58). 

A customary representation for the surface wave cross section functions 


is 


*£■" (y) 
*r en (y) 


M^w-r/^AS -) 1 ; |y| > r/2 



; |y| > r/2 


(61) 

(62) 
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where \p£ ven < odd must now be chosen so that they satisfy their associated 
GSTC. By substituting (61) and (62) into (10) and (11), we find that 
A odd.even mus t satisfy the polynomial equations 

\fl - (A !!”“)’«& ([A5““] J ) + ttf, ([A'-"] 2 ) =0 (63) 


v'i -(x';r m ru' n +«;■, ((vr"f ) = o (64) 

and can be also identified as the poles of the slab plane wa r e reflection 
coefficient. We further note that 




foP^ven 

^ m 


(y) 


tyP.odd ^ 

\p,odd 

\ y ) 

U\2 ( 

[A P rn° dd ' 

’) 


^p.even ([A^ en ] 2 , y 

U r „ ([AS'”"] 2 ) 


|y| > r / 2 


\y\ > T l- 


(65) 


( 66 ) 


implying that for a multilayer slab the cross section functions associated 
with the discrete and continuous eigenmodes are of the same generic form 
given by (59) and (60). 


Interior Cross Section Functions 


We consider now the determination of the cross section functions for the 
region interior to the slab (i.e. in the region \y\ < r/2). For simplicity let 
us first assume a single layer slab of thickness r = r x , whose upper face is 
located at y = — ti/ 2 . In accordance with the preceeding, the cross section 
functions associated with the external fields are given by 

(Wy) = M {?„(«;,<, rf, A 2 ) i, a (l,l,|y|-r./2, A 2 ) 

+ „ 2 (u;.*;,Tf,A 2 ) fe(l,l,ll/|-T./2,A*)}; |y| > t,/2 (67) 

*—{x\y) = {to(<*f,rf,A 2 )5 12 (l,l,l!/|-n/2,A 2 ) 

+?22 (“i. rf, A 2 ) 933(1, 1, |y| — Ti/ 2, A 2 )} ; jy| > r,/2 (68) 
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obtained by setting Ufj (A 2 ) = y,j (u p , k\, rf, A 2 ) in (59) and (60). These 
are orthogonal functions and each must, therefore, satisfy the continuity 
conditions 


y P ,odd (aVj-) 

—dy* p ’ odd (A 2 ,rr) 

Ui v ' 

^ Pt e,en (A 2 , r") 

—dyW' even (A 2 , r r) 

U\ V ' 

with similar conditions on odd, even/ 

m Vi 

that possible cross section functions 


= tf p ’ 0<w (a 2 , r+) 

(69) 

= dy<b p ’ odd (A 2 ,r, + ) 

(70) 

= 'P p ' even (A 2 ,r+) 

(71) 

= dy'I' p,even (a 2 , Tj + 

) (72) 

). It is now straightforward to deduce 
satisfying (69) - (72) are of the form 


* pMd (a 2 , y) = (73) 

'K'""" (A 2 , y) = »,(»?. A 2 ) (74) 

for |y| < r j 2. Also, in view of (69) - (72), the cross section functions for 
the surface wave modes remain as given in (65) - (66), provided (73) and 
(74) are used in place of \^ odd ^ ven # 

For the general case of a multilayer slab, it is necessary that each of 
the internal cross sections functions satisfy the continuity conditions at all 
layer interfaces comprising the slab. In this case we find that 


p i odd 




^fi(A 2 ) 912(1, l,|y| -t/ 2, A 2 ) 
+^ 2 (A 2 )g 22 (l,l,M-r/2,A 2 ); 

|y| > r/2 

Vn (A 2 )?i 2 (wf,«f,|y| -y/_i,A 2 ) 
+'P l i2 (A 2 ) ? 22 («f, #cf, |y| - y/_i , A 2 ); 

y/ > |y| > y/-i 
912(^1, |y|. A 2 ); |y| < yi 


(75) 
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yp.even ^2 ^ 


< 


(76) 


where 


^2 P i(A 2 )?i2(l,l,|y|-r/2,A 2 ) 

+7/22 ( A 2 ) <722(1, 1 , |y| — r/ 2 , A 2 ); 

|y| > r/2 

Pjf (A 2 )y 12 (uf,/cf, |y| - y/_i, A 2 ) 
+7 > 22 (A 2 ) ? 22 (wf, «f, |y| — y;-i, A 2 ); 

y/ > |y| > y;-i 

, 922(w?,«i, |y|,A 2 ); |y| < yi 


T)^P / 3j: 2 \ W.P 7 \ 

^21 (“ k 2 ) ^22 


=n 

m=l 


912(m 


<, r£ , -*-£) y 2 2«, t£, -^) 


*.-p 


t p _ 

? ' m i ^2 / 

Ox 2 


(77 


When these are used in (53) - (56) in conjunction with (65) and (66) we 
have a complete field representation for all x. 


6 Recasting of the Dual Integral Equation 
Solution for a Material Junction 

The expressions (53) - (56) can be used to represent the fields interior and 
exterior to the slab. It remains to find the coefficients of these expansions 
and to do this we must first rewrite F(x,y) in a form compatible with (53) 
- (56). That is, we need to identify from (13) and (14) the discrete and 
continuous spectral components. The discrete portion of the spectrum is, 
of course, comprised of the geometrical optics and the surface wave fields. 
These can be identified by detouring the integration path in (14) as shown 
in Figure 7. In particular, for x < 0 the integration path may be deformed 
to one over the branch cut in the upper half of the A plane, capturing 
any surface wave poles attributed to the zeros of ^^(A) and Q\ v * n (A). 
Similarly, for x > 0, the integration path may be deformed to one over 
the branch cut in the lower half of the A plane causing the capture of the 
geometrical optics pole at A = —\ 0 in addition to any surface wave poles 
attributed to the zeros of (A) and (A). 
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Through the above deformation of the integration paths in (14) we 
obtain 

P o M( , _ ( F}f“ (*, y) + Fl£ dd (*, y) + fij? (*, y) x < o 

1 ’ y) l (*, y) + ***" (x, y) + (. X , y) x > 0 ( j 


—ven, v = f Fj/ ven (x,y) + F^ eve "(x,y) + / j t -;7 n (x,y) x < 0 

1 Fgo ven (*, y ) + F»w ven (x, y) + Fdi} v f n (*» y) x > 0 


0) 


where the components F ff0 , F sw , F*-// denote the geometrical optics, surface 
wave, and branch cut (or diffraction) contributions to the total fields. 
After some manipulation we find 

tx COS 4>o 


Fgo° dd (x, y) = A\' odd (\ 0 )* hodd (\ly)e> kl 
F]f d (x, y) = Al odd (A 0 ) T 2 ' 0<w (A 2 , y) 

F' g ; even (x, y) = J 4 j’ et ’ en (a o ) ’J ,1,euen (A 2 ,y) e j kxcos<t, ° 
F 2 30 even (x, y) = -Ai’ even (A 0 ) ^ 2 ’ ewen (a 2 , y) e J kxcos ‘ t ‘° 

where the A expansion coefficients are identified as 

sin 4> 0 e jfcT / 2sin *° 


.4}’ 0<w (A o ) = 


.4 2l0 "(A 0 ) = 
Ap* en (A 0 ) = 
^’ ewen (A 0 ) = 

For the I? coefficients, we have 






1 ,odd 


E B, 0.) 

/=1 


Q f"(A|) 

sin^„e J,:T/2 ” n *» 

ar“(A*) 

sin <f> 0 e jkT / 2sia ‘t > ° 

Q{ ven {K) 

sin <j) 0 e 3kT l 2s ' n *° 

Qr n ( A 2 ) 




(SO) 

( 81 ) 

( 82 ) 

(83) 

(84) 

(85) 

( 86 ) 
(87) 
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( 88 ) 


N odd~ 1 ^odd — 1 - m 

S amn 

m = l n=0 


(x}' oJd 




^ 1 ,odd 



e 


— jkx\ 


l ,odd 
( 


^1 ,ei/en 


with 


B]' odd {\ 0 ) 


1 .even 


{i,y) 


\ B, A,) 

/=1 

^ e u « n “ 1 ^even“ l ” TO 


£ £ b„ n (x^ + k)’" (x'r” k)" 


m = 1 


n=0 


^l.efen ( [A|' et ' en |\y) e'^^' 


1 ,et/en 


(89) 


— sin <p 0 

a! ,d<w + A 0 


D jkr /2 sin <£ 0 


3 £,°£ d ( A ) 

al 


J A=.\ 


^(Ao)^t'(A/ 1,0 ") W(A 0 ) 


(90) 


B\' cvcn (A 0 ) = . Sm< ^° 

f V ° f \ 1 ,ei>en . \ 

■A/ “t" 


pjkr /2 sin <J > 0 


wrr n (A)i 

aA 


J A=A 


(9D 

af- n (Ao) ^2 e + " (A?’ euen ) £!+" ( A*) 


The expressions for F^ dd (x, y) and F^ ven (x,y) parallel those in (86) and 

(87) . 

To obtain the C coefficient we express Fj-Jj as 


777 1 ,odd 

t di! } 


{*,y) 


y_ f°° J_ -ffy / 1 ~ A^ 

\y\ 7-oo 2n yi _ ^2 _ ^2 _j_ A o j 


^-jk\y\0 g-jkxy/TFfF 

"ot d Jyf^WTs^) of- d (K)Q° 2 d + d (yr^) ^(Ao) 
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where the branch of the square root is chosen so that /■ y/1 — /i 2 ) > 0 and 

6 is a vanishingly small positive number. By splitting the integral into its 
positive and negative portions, and employing some identities (92) becomes 

= jH C hodd (0) Z odd (^/l — /? 2 , A 0 ) 

. (i _ /3 2 , (93) 


with the expansion coefficient C l ' odd (0) given by 



Similar expressions for F%-j dd (x, y ), Fj-j V j en (x, y ) and F^J 71 can be obtained 
in a parallel manner leading to the identification of the remaining C coef- 
ficient. 


7 Determination of the Constants 

To determine the constants a mn and b mn , we may now enforce the tangential 
field continuity conditions 

F(x = 0“,y) = = 0 + ,y) ; \y\ < r/2 (95) 

dxF(x,y ) x=0 . _ = -±-jdxF(x,y) x=0+ ; |y| < r/2 (96) 
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with 


«i,2 (y) 


A*i ,2 (y) E z - pol 
«u(v) ^-pol 


(97) 


and the subscripts 1 and 2 denoting quantities attributed to the left and 
right side of the slab. Substituting (79) - (79) into (95) and (96), we obtain 


Fgo° dd (x = O', y) + Flf d (x = 0 ”, y) + (x = 0", y) 

= F]f d (x = 0 + , y) + F,£ dd (x = 0 + , y) + (x = 0 + , y) (98) 


ui (y) 


ax [F^“(X.I,) + F! w ° M (x,y) + Fir// (*,!/)]„„_ 

= + F% M (x,y) + Flff (x,y)] 


x=0+ 


(99) 


(x = 0- , y) + fir- (* = 0-, y) + F];?r (x = 0- , y) 

= FjT"(x = 0 + ,y) + F'S’'{x = 0 + ,y) + F$T(x = 0+,y) (100) 


^ax (x,y) + F'S m (x, y) + Fi>7” (x,y)j i=0 _ 

= ^ax[^rUx,y) + F'r”(x,y) + F^7”(x,y)] i=0 , (101) 


to be solved for all a mn and b mn . In particular, for an odd GSTC of 0(N^ dd ) 
to the left and of 0(N% dd ) to the right of the discontinuity, the number of 
a mn to be determined is equal to 


N a = 


N 0 dd (NoJJ — l) 

2 

( ‘^dd + ^Idd ) ( ‘^Idd ~t ^odd ~ 2 ) . 

I 8 > 


N odd + N odd is even 

N] id + Nl dd is odd 


(102) 


To determine all a constants, (99) and/or (100) must then be enforced or 
sampled at a minimum of N a points accross |y| < r/2 and 0 < <f> 0 < tt. 
Similarly for an even GSTC of 0(N* ven ) to the left and of 0(jV| t ' en ) to the 
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right of the discontinuity, 


N b = 


^ even (i ^ even 


(^.n + ^«n)(^n + ^\«n-2) . pyl , ^2 

S i ^ even ' ^ ev\ 


is even 


( /V 1 +/V 2 ) -1 

^ 1 even ' 1 T even / 1 

8 


TV 1 4- N 2 is odd 

v even ' ^ v et;en uuu 


(103) 


and thus, the b constants can be determined by enforcing (101) and/or 
(101) at a minim. of N \ points. 

Substituting ior the fields in (99) and (100) as given in the previous 
section, we obtain the equations 


Vp dd (A 0 , y) = ^ a p Zp d (m (p ) , n (p) , A 0 , y) 

P-1 

VftS(A..lf) = £«p-Z s 1f (m(p),n(p),A„!,) 

P-1 


where a p = a m(p)n(p ) with 


(n + m — l) (m + n) 

p = +m 


, x P ,f \/l +8(p - 1) - 1 

m {p) = p — —Int 1 


• Int 


n(p) = Int 


f ^/l +8(p- 1) + 1 j 

yi+8( P -i)+i 


- m (p) 


(104) 

(105) 


(106) 


(107) 


(108) 


which are in accordance with the ordering of the a mn constants as the order 
of the GSTC is increased (see Figure 8). The functions V F dd , V£ d p, Z F dd and 
Zg dd F are readily determined from the previous analysis and are not quoted 
here. 

Equations similar to (104) - (105) can be obtained for the b constants 
in a parallel manner. 
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8 Validation of the Solution 


The validity of the derived angular spectra and diffraction coefficients was 
already performed to a limited degree in Section 4 of the paper. What 
remains, therefore, is a validation of the procedure for computing the con- 
stants a mn and b mn which amounts to solving a small matrix. The valida- 
tion was done by comparison with processed data from a numerical model 
which consisted of a finite length slab having the prescribed discontinuity 
or junction at its center. First, the transient response of this finite slab 
was generated from bandlimited frequency domain data. The contribution 
from the material junction was then obtained by time gating the transient 
response. Numerically derived data from this procedure were found in good 
agreement with the presented analytical solution. An example is shown in 
Figure 9 corresponding to a thick (0.2 freespace wavelengths) material half 
plane. The numerical and analytical data are clearly in good agreement, 
and it is again demonstrated that the constants play a major role in the 
solution. 

9 Other Applications of the GIBC/GSTC 

It was shown above that the GIBC can effectively model thick planar layers 
of material. However, corresponding GIBC can also be derived for curved 
coated surfaces (see Figure 10), and the improved accuracy of these is 
particularly evident when surface wave effects are dominant. For surfaces 
having relatively large radii of curvature these can be easily derived from 
those of the planar surface with x and z replaced by the local tangential 
variables and y by the normal one. With a second order GIBC derived in 
this manner, the Mie series and GTD solutions have been found [20], [21] 
and compared with the exact modal series solution for a coated cylinder. 
As illustrated in Figure 11, the field given by the Mie series based on the 
GIBC is in excellent agreement with the exact result even at points close to 
the surface in the shadow region where a finite order boundary condition is 
inadequate. In contrast, data based on the standard impedance boundary 
condition (SIBC) are substantially inaccurate. 

Higher order boundary conditions have advantages in numerical treat- 
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ments as well. When used to simulate a coating, a GIBC eliminates the 
need to sample inside the dielectric, and this is important when storage is 
limited. In addition, it may be possible to use a GIBC to transfer a bound- 
ary condition to a plane, thereby producing a boundary integral equation of 
convolution type. In conjunction with an FFT, the equation can be solved 
iteratively to reduce the storage requirement to 0(n) where n is the num- 
ber of unknowns. As an example, for the three dimensional problem of a 
cavity in a coated ground plane, a GIBC provides a simple modal as well as 
a reduction in memory. If the coatm^ ;s lossy or tapered in thickness, the 
non-uniqueness due to the terminations is avoided, and the same is true for 
cavities whose depth tapers to zero. Nevertheless, caution must be exer- 
cised when solving the integral equation numerically. The GIBC results in 
higher order derivatives applied to the Green’s function, and even if some 
can be transferred to the current, the increased singularity of the Green’s 
function makes discretization more difficult. In spite of this, integral equa- 
tion methods using GIBCs of up to the third oder have been successfully 
implemented [22]. 
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Appendix: Multiplicative Split Functions 

In this appendix we consider the splitting of 

Q (a 2 ) = U A (A 2 ) + v/T - \ 2 U b (a 2 ) (109) 

as a product of two functions, one of which is free of poles, zeros and 
branch cuts in the upper half of the A plane and the other having the same 
properties in the lower half of the A plane. That is, we seek to write Q (A 2 ) 
in the form 

C7 (A 2 ) = (A)^_ (A) (110) 

where the superscript + and - indicate an upper or lower function, respec- 
tively. Noting that 


U A (a 2 ) 

= E-Ml-Af 

(in 

n=0 


Ub (A 2 ) 

Nb 

n=0 

(112 


with N \ = Nb or Na = N& + 1, we may rewrite Q (A 2 ) as 

N s 

a (A 2 ) = g5 n [v^] n (113) 

n=0 

where N 3 = Mcix(2Na,2Nb + 1 ) and S n ~ A n j 2 if n is even and S n = 
B ( n -\)/2 if n is odd. However, since we seek a multiplicative splitting of 
(113), a more convenient form to represent Q (A 2 ) is 

g (a 2 ) = So n ( : + (H4) 

in which 7 n denote the zeros of the polynomial Ylhlo We immedi- 

ately now identify that each of the product terms in (114) can be factored 
as 

1 + vTET d g M (A; 7) M- (A; 7) (115) 

7 
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where 


(116) 


A'+(A;l/7) 


y/1 — A 
M+( A; 7 ) 


is the split function characteristic to the impedance half plane having a 
constant surface impedance I /7 [17]. With the branch choosen so that 
/m(V 1 — A 2 ) < 0 , M+ (A; 7 ) is explicitly given by 


M+ (A; 7 ) = M- ( — A; 7 ) = 


M + (A; 7 ) 

jn 

A/ + (A;--y) 


/m( 7 ) < 0 
I n(y) > 0 , 


(117) 


M+ ( cos a; 1 / ) 


n (tt/ 2) [l + y/2cos (^^)] [l + a/ 2 cos (ggg=^. 

/J [*„ (3tt/ 2 - a - 0) , (tt/ 2 - a + 0)} 2 

(118) 


In this, 


Im(rj) > 
A = 

Im (y]\ - 1 /t 7 2 ) < 

0 = 


0 

cos a: 

0 

sin _ 1 ( 7 ;) with 0 < Re(0), 


( 119 ) 


and 'h*. (a) is the Maliuzhinets function [18] whose evaluation in algebraic 
form has been given in [19]. 

The determination of Q± (A) is now rather trivial. By substituting (115) 
into (114) we easily obtain 

J,(») = !!.(-») = l/s>n«t(UJ (120) 

71=1 
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Figure 1. (a) Distributed sheet, (b) Distributed sheet discontinuity. 



Figure 2. Geometry of the material-material join. 



Figure 3. Illustration of the C contour in the complex X-plane. 




Figure 5. Recessed slab (PEC stub) on a ground plane. 
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Figure 6. H z -polarization backscatter echowidth for a material insen having w=0.04X, 
e=2-j.0001, ji= 1.2 modeled with O(w) low contrast GIBCs. 
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Figure 7. Deformation of the C contour for (a) region l integrals and (b) region 2 
integrals. 
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Figure 8. Indexing scheme for the constants. 
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Figure 9. H z -polarization backscatter echowidth for material half-plane with t=.20, 
e=2-j.0001, |i=1.2. 





Figure 10. Illustration of a three-layer coated cylinder. 
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Figure 1 1 . Bistatic H-polarization scattering pattern of a circular cylinder of Radius 2.931 
coated with a layer 0.07A. thick having e=4 and fl=l ; Comparison of fields at a distance 
0.05a, away from the coatings surface. 





Table 2: Even symmetry parameters associated with the right hand side material 
Appendix for the definition of the A/+ split functions. 



















